Quantum teleportation of nonclassical wave-packets: An effective multimode theory 
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We develop a simple and efficient theoretical model to understand the quantum properties of broadband con- 
tinuous variable quantum teleportation. We show that, if stated properly, the problem of multimode teleportation 
can be simplified to teleportation of a single effective mode that describes the input state temporal characteristic. 
Using that model, we show how the finite bandwidth of squeezing and external noise in the classical channel 
affect the output teleported quantum field. We choose an approach that is especially relevant for the case of 
non-Gaussian non-classical quantum states and we finally back-test our model with recent experimental results. 
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I. INTRODUCTION 

As a striking example of quantum communication protocol, 
teleportation was discovered early on in the development of 
the field of quantum information processing. With either qubit 
ilt] and continuous variable flavors [2], experiments were soon 
to follow [3ii 4J- Until now continuous variable teleportation 
has only been performed with the class of so-called Gaussian 
states [4i-6]. However, this alone is not sufficient for universal 
quantum computation where non-Gaussianity of some kind 
has been shown to be necessary [7]. Although non-Gaussian 
nonclassical states of light that would allow for such universal 
operations have been available experimentally for some time 
in the continuous variable regime |8- 10], the major challenge 
of actually manipulating these states in some Gaussian pro- 
tocol context beyond simple generation has remained mostly 
unaddressed. 

Some recent experimental work has reported on success- 
ful continuous variable teleportation of a strongly nonclassi- 
cal state of light ifTlll . In this experiment, a close approxima- 
tion of a Schroedinger's cat state generated with the photon- 
subtraction protocol [12:1 is sent through a continuous variable 
teleporter. The quality of teleportation is high enough that the 
output teleported state is also a nonclassical state with a nega- 
tive Wigner function. The use of a non-Gaussian nonclassical 
state as an input state is the first most noticeable feature of this 
experiment. Although Gaussian states teleportation has been 
amply studied, due to the complex nature of non-Gaussian 
states and especially mixed non-Gaussian states, only few 
general results exist for this case. Maybe the most general 
condition for successful teleportation of non-Gaussian non- 
classical states is the necessary but not sufficient 2/3 thresh- 
old on fidelity |13]. In [141 . the theoretical work closest to 
the experimental conditions of iflTll . teleportation success is 
investigated for the case of a mixture of vacuum |0) and one 
photon |1) as an input state. On top of these difficulties, to 
accommodate with the transient nature of the input state used, 
the teleporter used in [11] operates on a broad range of fre- 
quencies. This is the second most noticeable feature of this 
experiment in contrast with typical continuous variable exper- 
iments, which only manipulate narrow sidebands of light. To 
our knowledge, there are actually very few results relevant to 
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the case of multimode teleportation. In 111511 . multimode tele- 
portation of a quantum field is investigated from the point of 
view of temporal fluctuations using the photon correlations 
function g^^^r). In 1 16], it is shown how to adapt the single- 
mode teleportation protocol of [ 17J to telepoitation of a multi- 
mode field with finite teleportation bandwidth. Both these 
works attack the problem of multi-mode teleportation from 
the Henseinberg picture and additional considerations are re- 
quired to handle the case of non-Gaussian nonclassical input 
states. 

Our main objective in this paper is to attempt to answer both 
these issues with a theoretical model as simple and efficient as 
possible. In Sec. [Ill we first briefly introduce the teleportation 
protocol with its usual phase-space formulation and define a 
criterion of success. In Sec. |III] we describe a realistic model 
of a nonclassical non-Gaussian state that faithfully models the 
input states of iflTll . With this model, we are able to predict the 
success of teleportation in a way similar to ]14]. In Sec. |IV] 
we use the Heisenberg picture to approach multimode telepor- 
tation as teleportation of a quantum field. We then show how 
to reduce this quantum field to a single effective mode that de- 
scribes the temporal properties of the input state. In Sec. |V] 
we show how to take into account any external noise spectrum 
in the broadband teleportation operation. Finally in Sec. |Vl] 
we compare our model with the recent experimental results of 
and conclude. 



II. BASICS 

Deciding on success of continuous variable teleportation is 
a non-trivial problem as it is closely related to the kind of 
input states and entanglement used, as well as the specific 
protocol or quantum circuit teleportation is actually used for. 
For the Gaussian case, the fidelity F = {ijjin\pout\'4'in) is 
the usual figure of merit, though F loses much of its mean- 
ing as a benchmark figure when more general non-Gaussian 
mixed states are used as input states. While Gaussian states 
can be fully characterized by their first and second moments, 
which allow figures like fidelity to have some general and 
useful meaning, such an approach fails with non-Gaussian 
states. Because non-classicality itself is an ambiguous prop- 
erty for continuous variable systems with infinite-dimensional 
Hilbert spaces, it is even more complex to decide on a relevant 
success criterion for continuous variable teleportation of non- 
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classical quantum states. In this paper we consider the input 
M^'" and output W°"^ Wigner functions of a teleportation pro- 
cess. We adopt as a criterion of success the successful transfer 
of negative features of the Wigner function. Provided is 
itself a negative Wigner function and having for W" a pre- 
cise algebraic expression including the relevant experimental 
parameters, we want to know what are the requirements on 
these parameters and on the teleportation process for success- 
ful retrieval of negativity in M^""'. Furthermore, we restrict 
ourselves to the Braunstein-Kimble scheme described in IItIi 
where the teleportation can be expressed in phase space as the 
following convolution 

W°"' = W'''oG,-., (2.1) 

with r the Einstein-Podolsky-Rosen (EPR) correlation param- 
eter and Ga {q, p) a normalized Gaussian of standard deviation 
a (h = 1). In this case, teleportation of nonclassical features 
such as negativity has been shown to require 3 dB of squeez- 
ing [Ts], or equivalently a vacuum fidelity of F > 2/3, which 
is also called the no-cloning limit [18]. Precisely speaking, 3- 
dB is a lower bound for unity gain teleportation of negativity 
of any pure or mixed state. Recent work has shown that, given 
the precise shape of the input state and amount of anticorre- 
lation in the teleportation quantum channel, there actually ex- 
ists strategies to surp ass the 3 dB threshold by tuning the gain 
of teleportation Ill4ll . However, unless extremely pure entan- 
glement is used, typical experimental antisqueezing imposes 
virtually unity gain operation. Furthermore, as this tuning is 
input dependent, the teleportation setup losses its universal 
characteristic. 

The 3-dB threshold is only a lower bound to negativity 
teleportation and we would like to have a model that pre- 
dicts better the success or failure of negativity teleportation. 
In the general case, this is too broad a problem to handle 
as teleportation is known to be input dependent. From now 
on, we focus our analysis on the specific case of the photon- 
subtracted squeezed vacuum that was used as an input state 
in [|1 1|1 . This family of quantum states has recently attracted 
a lot of interest both experimentally lIsl-IToll and theoretically 
lfl2i[T9ll20ll . From an optical parametric oscillator (OPO) pro- 
ducing a squeezed vacuum Ss\0) with squeezing parameter s 
called the signal mode, a fraction R of the output called the 
trigger mode is tapped and sent to a photon resolving detec- 
tor to herald non-Gaussian states. Various single-m ode and 
multi-mode models exist for this protocol lll2l Il9l I20I1 . and 
they are essentially equivalent in the limit of small s and R. 
To start, we will assume that a detection event projects Ss\0) 
on the photon-subtracted squeezed vacuum aiSsjO) equal to 
a squeezed photon 5*8 The Wigner function W'^'^^ of this 
reference state is written 

W'^\q,p) = 2{e^^q^ + e-V - l/2)Gi/^(e^g, fi-^p), 

(2.2) 

and has a maximal central negativity of W''^^{Q, 0) = — I/tt. 
Although this is a specificity of this particular input state, from 
now on we will use the value of the Wigner function at the 
origin of phase space as the figure of merit for negativity tele- 
portation. Applying teleportation equation (12.1b on Eq. (12.21 ). 



we find for output negativity Wref out(0, 0) 

7r((l + 2e-2'-)+8e-2'-sh2(s))-'^ 

which indeed yields W,ef out(0, 0) < for r > In -\/2 as ex- 
pected. IVref out(0, 0) will become negative only if the param- 
eter r is greater than In \/2, equivalent to 3 dB of squeezing. 

III. REALISTIC INPUT STATE 

Equation ( 12.31 1 is, of course, of little interest since an actual 
experimental input state will virtually be a mixed state and 
a more realistic model of input is required. The experimen- 
tal input state used in |11] happens to fit well a simple loss 
model where the experimental Wigner function M^'" can be 
modeled from W^^, the reference state, by applying "beam- 
splitter losses" 1— 77 equivalent in phase space to the operation: 

W'\x,p) = i (M^-f o Ga) (^-^, ^ , (3.1) 

with A = \/^^r^ |21]. The phase space transformation (13. Il l 

is derived from the action of a fictitious beam splitter with 
transmission coefficient 77, which transforms input coherent 
states I a) into I77Q;). As was shown in |22], it is also possible to 
express the action of this beamsplitter with a master equation 
acting on the density matrix p. In this case, we would obtain 

^p(t) = L[p{t)], L[p] = [pa,a^ - [pa\a]. (3.2) 

By using the previous algebraic expression (12.2) of W'^^ 
together with the transformation (I3.lt . we can obtain the exact 
expression of W™ with any mathematical software. We rather 
are interested in the central negativity given by 

iy™(0, 0) = (1 - 277)/^ (1 + 477(1 - 7/)sh2(s))^/' , (3.3) 

where the negativity threshold M^'"(0, 0) — depends only 
on 77: 77 > 0.5 implies W™{Q,0) < (see Fig. [T]). Using 
this model of imperfect input state we investigate negativity 
teleportation of W" by concatenating Eqs. ( 13.11 ) and ( 12.11 ). 
The two successive Gaussian convolutions are reduced to one, 
while the phase space rescaling survives the teleportation. 
Eventually 14^°"' happens to be written in the same form as 




(3.4) 



where A has been changed to A' = y A^ + (e""")^ /rj in a 
way similar to classical amplifiers input/output SNR rules. 
We remark that besides the degradation of the input nega- 
tivity VF'"(0,0), the loss parameter 1 — 77 has also the ef- 
fect of decreasing the effective correlation parameter r to 
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Figure 1. Input Wigner function negativity W"(0, 0) as a function 
of rj for different values of tfie squeezing parameter s. 
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Figure 2. Output Wigner function negativity W^'"(0, 0) as a function 
of r for different values of rj and s. For a given rj, different s curves 
cross the W{0, 0) = at the same r{ri). 



r' = r + In ^ < r. In practical terms, this means that both 
operations do not commute and losses at the input stage have 
more effect on the quality of the overall process than losses 
at the output stage. The output center negativity is now ex- 
pressed as 



T^°"'(0,0) 



grigr - 2??) 



3/2 ' 



(3.5) 



with gr = 1 + 2e~^'' (plotted in Fig. |2]i. As expected for unity 
gain teleportation, the iy°'"(0, 0) = threshold is still inde- 
pendent of the squeezing parameter s and can be expressed as 
a function of the two parameters i] and r alone by the simple 
relation (Fig. [3j 



r In \/2/{2ri - 1) at threshold. 



(3.6) 
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Figure 3. Negativity threshold as a function of r and rj. 

State with an initial density matrix p^^f of the form 
A-ef = a5s|0)(0|5_sa'''. 



(3.7) 



However, to be more faithful to experimentally produced 
squeezed photon states, rather than an ideal photon resolving 
detector, we should consider the unideal projection properties 
of the Geiger silicon Avalanche PhotoDiode (APDs) experi- 
mentally used to produce photon subtracted states. Essentially 
two mixing mechanisms are at work. First, the on/off charac- 
ter of the APD makes it only able to detect the presence of 
photons without resolving the actual number of them. This 
leads quite naturally to a Positive Operator Valued Measure 
(POVM) solution to model the APD measurement, as was 
done in ll23ll . However, the effect is rather marginal if we 
restrict ourselves to small squeezing s and small tapping frac- 
tion R and we will actually neglect photon components higher 
than n = 1 in the trigger channel. Second, the laboratory APD 
is also characterized by a dark count rate, which will produce 
false heralding events and induce some statistical mixing of 
the target state. When such a false event occurs, no projection 
happens on the signal mode as no photon is subtracted and the 
signal mode density matrix is just pfaise = S's|0)(0|S'-s corre- 
sponding to the OPO output squeezed vacuum. We introduce 
a parameter e that reflects this statistical mixing and write the 
new reference density matrix p[.^f as 



Pref = (1 ^ £)A-ef + C/Ofalse, 



(3.8) 



where e is related to the modal purity parameter 5 introduced 
in [24] by 1 - e = S. With e ^ 0, the value of W^'"'=f(0, 0) is 
not optimal anymore, but becomes (Fig. |4|i 



W^"=*(0, 0) ^(1 - e).VK'"'=*(0, 0) + e.H/'^'^'=(0, 0) 
-(2e-l)/7r. 



(3.9) 



Correcting for the effect of Eq. (13.8) in the input negativity 
(|3.3l l, output negativity ( 13.51 ), and negativity threshold (|3.6l l is 
just a matter of calculating how the Wigner function associ- 
ated to pfaise, WMse{x,p) — G^/^(e^x, 6""^), cvolvcs in the 
teleportation process. This is possible since all the transfor- 
mations used until now have been linear and therefore we can 
write 



M^"(0,0) -> (1 - e)W^"(0,0) 



Until now, our model of reference state has assumed a pure 



+ -(malseOGA»)(0,0) 
77 



(3.10) 
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Figure 4. (a) Input negativity threshold V(^'°(0, 0) = as a function 
of Tj, s, and e. (b) Input negativity M^'"(0, 0) as a function of rj for 
different values of s and e. 



with a e {in, out}. A'" = A, and A°"' = A'. We find for the 
corrected input negativity the new expression 

T^'"(0,0) ^ M/'"(0,0) 

l + 2(l-?7)sh^(s) 



7r(l + 4r;(l-7?)sh2(s))^^^' 



(3.11) 



as well as 



W^™'(0,0) ^ W^™'(0,0) 

,0 9r + 2{gr - T?)sh^(5) 
+2ejj ryj, (3.12) 

71(52 + 4,7(5, -2r^)sh2(s))'/' 

for the corrected output negativity. Input state negativity 
threshold VK'"(0, 0) = now gives the following relation be- 
tween T] and e 



= (2,7-l)/277(l + 2(l-7?)sh2(,)) 



(3.13) 



For nonzero e, the input threshold becomes dependent on the 
squeezing parameter s. We expect the same dependence on 
the output negativity threshold corrected for e, which is now 
expressed with the following quadratic equation: 



= gl + 2b{e,s)Tjgr-cie,s)7]^ 



(3.14) 



with 6(e,s) = e(l + 2sh^(s)) - 1 and c(e, s ) = 4esh^(s). 
Keeping the only physical solution in ( 13.141 ) the negativity 
teleportation threshold corrected for e becomes (Fig. |5]l 



1/2 



r = In 



Tj{^/WT7:-b) - 1^ 



(3.15) 



In summary, we have developed in this section a realistic 
yet simple model to account for the phase-space properties of 
photon subtracted squeezed vacuum states. With that model 
we have considered the effect of teleportation on these states 
and are able to predict the success of negativity teleportation. 
In the next section we will now show how to take into account 
the multi-mode aspect of these input states and the multi-mode 
aspect of broadband teleportation. 
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Figure 5. Influence of e and s on the negativity tftreshold 
W°"\0, 0) = at the output of teleportation. 



IV. MULTIMODE TELEPORTATION 

The APD triggered non-Gaussian state has been shown to 
have complex multimode properties ||20[ I25I1 . On one hand, 
the OPO output beam is a continuous wave with a specific 
squeezing spectrum, while, on the other hand, the APD trig- 
gers that herald a non-Gaussian state happen at precisely de- 
fined times. Intuitively one expects that a given APD trigger 
induces non-Gaussian statistics in the immediate time vicinity 
of the trigger event to the extent of the OPO bandwidth. When 
the number of triggering events increases as R and s increase 
and more photons end up in the trigger channel, more com- 
plex time interference phenomena arise on the signal mode be- 
tween neighboring photon subtracted wave packets 1 26 J. For- 
tunately, in the limit of small s and small R, a simple two 
modes picture allows one to efficiently describe the input state 
and capture most of its experimental properties [19J. It in- 
volves two effective wave-packet modes, Ag for the signal 
mode and At for the trigger mode, defined by 



A, 



(4.1) 



with i G {s, t}. Preserving the commutators [Ai, A]] = 1 re- 
quires J \fi{uj)\'^duj = 1. While the exact form of ft{oj) is not 
really relevant, since the APD detection time is typically much 
shorter than any other time scale in these continuous wave ex- 
periments, /s(w) will describe the temporal characteristic of 
the heralded non-Gaussian state. This function will be defined 
by the OPO bandwidth, as well as the possible filtering cav- 
ities used on the trigger channel and numerical optimization 
has shown that the optimal form (small s and R, wide filter- 
ing cavities) can be taken as 



/(cj) = 7/7r[7^ + {ojQ - w)^ 



(4.2) 



with 7 the OPO decay rate including^ntra-cavity losses and 
LjQ the light beam carrier frequency[26|. In short, an APD 
heralded state behaves as a traveling wave packet of light with 
non-Gaussian characteristics. Using the definition of As in 
Eq. ( 14.1b . we return to a pure single-mode model for the input 
state that we write as 



-s{A^' 



')/2it|0). 



(4.3) 
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This expression is still an approximation of reality in the sense 
that a true multimode description as done in |20J would re- 
quire an extension of a broadband squeezing operator 5'^ on a 
basis of adequately chosen orthogonal functions {/„(tj)}. In 
this case, Sb is expressed as 



Sb = exp 



2^ 2 



AnA- 



(4.4) 



with Afi = a{ujQ + fi). In practice, once a signal mode 
fs = fo is chosen, the other modes n 7^ are traced out, 
which leads to mixing of the density matrix to the extent of 
the multimode entanglement present in Sb between all modes 



Pb = tr„5^o 



(4.5) 



This multimode entanglement is quantitatively tracked by the 
function C{^)- A special case happens when is a constant at 
every frequency. Then the operator Sb can be exactly factor- 
ized on the orthogonal basis {/„(cj)} and the expression i4.3i 
becomes exact. In the general case of a non constant func- 
tion the expression (14.3b is nevertheless useful as it only 
neglects a small amount of entanglement between the differ- 
ent orthogonal modes n if s is small. As a result, we will 
use expressions ( 14.1b and (14.3b for our following analysis of 
multimode teleportation. 

To investigate how this multimode aspect translates quan- 
titatively, we express in the Heisenberg picture the relation 
between input {xin,Pin) and output (xouuPout) quadrature op- 
erators of teleportation as 



-e H T^e^ Wx, (4.6) 



V2 



Pout = gpPin H ^1^^ 



V2 

^e+'^f),, (4.7) 



with (vx , Vp) and {lix , lip) two auxiliary modes in the vacuum 
state 1 27]. We first consider the unity gain case gx ~ gp ~ 1, 
where the input/output relations (14. 6b and (14. 7b simplify to 



V2( 



e Vx , Pout = Pin + 



V2e ''wr, 



(4.8) 



We notice that the input modes Xi^, Pin and output modes Xout, 
Pout can actually represent any frequency mode cj and we de- 
fine in the same way as in Eq. (14.11 ) two new modes, the input 
Ain and output Aout wave-packet modes by 



Ain = J /s(w)ain(a;)dw, Aout = J /s(a;)aout(t^)rfw. (4.9) 

as well as the input and output wave-packet quadratures 
{XinjPin) and (Xouti^out) relevant for wave-packet telepor- 
tation. We can directly rewrite the teleportation input/output 
relationship ( 14.8b in the form 

^out ^Xin-V2j /,(c^)e-''(")fi, duj, (4.10) 
-Pout ^nn + V2j f,iu;)e'^'-'^^Wpduj, (4.11) 



where we have introduced r(a;) as the spectrum of EPR cor- 
relations resolved in frequency. From the physical properties 
of the OPO cavities used for EPR squeezing generation it is 
possible to deduce the expression of r{uj) from the squeezing 

spectrum = (A'^xl^Auj)A'^x^„d(i^))/(^'^x^ac)- See, 

for example, Ref. lUol for details. We now define an effective 
broadband EPR parameter r^ff by 

e-"-'"^ //,(L^)e-'-(")dw= f fsiLu)S.{uj)duj, (4.12) 



and therefore Eqs. (14.10b simplify themselves to 

^out - ^in - V2e-'-«» Vx, Pout - Pn, + V2e~'^"' Wp. (4.13) 

Since the auxiliary modes Vx and Wp are effectively traced 
out on the vacuum state at all frequencies, it is possible to 
take them outside of the frequency domain integrals to obtain 
the formulation (14.13b . We notice that Eqs. (14.13b are writ- 
ten in the same way as Eqs. (14.81 ). Thanks to the linearity of 
transformation (14. 8b and the linear model of input state (14. 11 ) 
in the Heisenberg picture, multimode teleportation is equiva- 
lent to familiar single-mode teleportation, where an effective 
broadband EPR parameter rgff has been defined to take into 
account the finite bandwidth of entanglement. In short, all the 
previous formulas of Sec. |III]for unity-gain teleportation are 
readily usable with the simple change r — > reff . 

The case of nonunity gain multimode teleportation is much 
more complex and we conclude this section with a brief 
overview of the non-unity gain case. First we introduce the 
transfer functions gx{ui) and gp{uj), which represent the ef- 
fects in frequency space of the classical channel, gx and gp are 
in general complex- valued functions verifying the Kramers - 
Kronig relations. We are now facing the problem that the 
output quadratures operators of teleportation will not be Her- 
mitian operators anymore, in general. By taking an approach 
similar to Eq. ( 14. 10b . we obtain for the position quadrature 



^o"ut = J fsiu^)9xic^)xU^)duj 

--^ (^e-^'"' + J fs{u;)gx{uj)e-^^"Uuj^ Vx 

+ (^e+''- - j fs{uj)gx{uj)e+^^"'>du?j Wx- (4.14) 

where is the output teleported mode in the nonunity 

gain regime. To clean this expression we define the two com- 
plex numbers : 



gt = I /.(c.).g,(u;)e±('^('^)-'--)dw, (4.15) 



so that Eq. ( 14. 14b simplifies to 



x^nonunit 
-^out 



fsiu!)gx{^)xin{uj)du: 



V2 



V2 



Wx, 



(4.16) 
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with a similar expression for P"°™"". By further separat- 
ing X" """"" and -P"°™"" in real and imaginary parts as done 
in il6ll . it is possible to obtain from this model observable 
results. Equations ( |4.16t and ( |4.15t show that the output 
modes will get contaminated by antisqueezing when [^^.(ti;)! 
and |.gp(w)| are different from 1. Furthermore, the expression 
/ fs {^)9x ('^)a;in {i^)duj hints that the wave-packet shape from 
input to output will get modified by the teleportation process. 
An interesting and practical situation is the case of pure linear 
delay gx{(^) = = exp[— iw Ai]. If such a phase factor 

is added by the classical channel to the output modes, its full 
effect can be absorbed in x-m and auxiliary modes Vx, Wx by 
using their Fourier transforms. 



-iuj At 



1 



dtXm(t) 



SO that is related to Xout by a simple time translation 



(4.17) 



As a matter of fact, this is exactly how the experimental tele- 
portation setup used in ifTTll behaves, where an optical delay 
line is used to match the phase answer of the classical channel 
and cancel this At phase factor. 

In this section, we have developed an efficient model of 
unity-gain multimode teleportation with the added benefit of 
being able to use all the results of the previous section. In the 
next section we will further use this model to investigate the 
effects of classical sources of noise in the classical channel 
and their effect on the teleportation process. 



V. NOISE MODEL 

In this section, we try to better understand the effect of clas- 
sical sources of noise on the teleportation process. This is an 
important point to consider as, compared to the single side- 
band regime, it is much harder to experimentally insulate from 
external noise a broad range of frequencies at the same time. 
To throw light on that issue, we first look for a master equation 
describing the effect of teleportation on the density matrix p. 
For that, we start by fully detailing Eq. (12. Il l with input W 
and output W': 



and the second derivative of W with respect to position x 



1 



1 f f dx'dp' 



{x,p) ^2 ''{x.p) a'^ J J 27rcr2 
x{x-x'yW(^^,y)e ^ 



(5.3) 



We immediately find the following differential equation for 

|w/(,,p;t) = ^(5*^'W)AW/^,,^,), (5.4) 

where A — + dp. We choose a to be a{t) — \j2K't with 
k! a constant decay rate, so that Eq. (15. 4t simplifies itself to a 
pure diffusion equation. 



n!AWl 



{x,p;t)- 



(5.5) 



Then, by using correspondence rules between the phase space 
formalism and the density matrix formalism ll28l . we find from 
Eq. i5.5i the following master equation for p: 



d_ 
dt' 



-p — k! {2a} pd + 2dpd^ 
—a) dp — dd^ p — pd^ d — pdd^^ , (5.6) 
which can be equivalently written as 

J^P = L[PI L[p] = K'[d\[p,d]]+K'[d,[p,d% (5.7) 

The master equation (15.6b is the well-known damping pro- 
cess for the harmonic oscillator Equations (14. 8t look there- 
fore similar to quantum Langevin equations, where the terms 
^/2e~'^Vx and ^/2e^'^Vp are nothing other than thermalization 
terms. We would intuitively add the effect of any classical 
source of noise directly in Eqs. (14. 8t by writing 



Pout 



= fin - V2e '"uj; 
= Pin + V2e^^Wp 



V2N'xy, 
- V27Vp z, 



(5.8) 



1 



dx'dp'W(x' ,p')e 



where we have introduced two new auxiliary vacuum modes 
y and z and where we first are considering the single-mode 
case. Afx and Afp describe the amplitude of noise normal- 
ized to vacuum and added at the output of teleportation on 
top of finite squeezing. This noise can arise independently 
for both quadratures from imperfect electronics in the classi- 
(^„^')2^.(p_p<)2 cal channel, for example. The most natural case is for noise 

■ to be uncorrected with quadrature angle and we can assume 

(5.1) -^x = A/'p = A/" to be the average noise amplitude. We re- 



We now assume that a — > a{t) has a time dependence. We 
express the first derivative of W' with respect to time t, 



^W' = — - 

dt '-^■P^ a{t) dt"" 

1 /■ /■ dx'dp 



member that all auxiliary modes appearing in the Heisenberg 
picture teleportation equations are traced out on the vacuum 
state and are uncorrected. Therefore it would be natural to 
redefine a correlation parameter r' modified by the amount of 
noise with the simple relation 



20-2 / / 27rcr2 



^(x',p') 



(5.9) 



x((x-x')' + b-y)')e- 



,(5.2) 



so that Eq. ( |5.8l l would be written as Eq. ( |4.8l l. However, this 
approach is wrong and the correlation parameter r' cannot be 
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redefined in amplitude but should be redefined in power by 
writing 

g-2r _^ g-2r' ^g-2r^_^2^ (5.10) 

It is possible to justify this expression rigorously by establish- 
ing the link between Heisenberg picture equations (14.8) and 
the original phase space formulation of Eq. ( 12. It . For that pur- 
pose, we introduce the characteristic function "Xia) related to 
the density matrix p by the Weyl expansion formula 

Xia) - tr(p^„) = (D^) , J daX{a)D-c.. (5.11) 

where Da is the displacement operator exp[aa^^ — a*a]. If we 
write a = {u + iv)/ \/2, then X ™d W are related by the 
following Fourier transform: 

W{x,p)^^J J dudvX{u,v)e''''-'-'''P. (5.12) 

Now if we consider the unitary transformation 

X x' — X — \/2^Vx, (5.13) 

with an auxiliary mode {vx,Vp) having commutators [x, v.j.] ~ 
[p,Vx] — and [wajjWp] = i, then the displacement operator 
Da is changed to 

Da-^ D'a^ Da(E) ^ Da <S) e+'^"'^*^ (5.14) 

with D" a displacement operator acting on mode v and a' = 
{—y/2j X u + i X 0)/\/2. To express the new characteristic 
function X'^ we have to evaluate the trace of D'" taken on the 
vacuum for mode v: 

tr {\0){0\DI,) = {0\DI,\0) = (5.15) 
This lead to the expression of X'i'^^ 

X'{u,v) = Xiu,v)e-''"^"/^ (5.16) 

which immediately translates to a Gaussian convolution, such 
as Eq. ( 12. Il l, for the Wigner function W having the Fourier 
relationship ( |5.12t between X ^nd W. In this case, we obtain 
the semiconvolution 

W'{x,p) = -L^ I da;'Ty(x',p)e-(=""^')/27^ (5.17) 
V27r7 J 

By also adding the transformation p ^ p' = p + V^jiip, we 
would finally obtain Eq. ( 12. Il l provided we define 7 as equal 
to cxp[— r]. If we now also consider the added noise term 
\f2Myn in Eqs. ( I5.8l l, X^iu, v) would be written 

X'{u,v)=X{u,v)e-'^'^'+^'^-'/\ (5.18) 

which justifies to redefine the correlation parameter r in power 
and not in amplitude as 

r ^ r' = r - In \/lTJPe^ . (5.19) 




in vacuum amplitude unit 
Figure 6. Influence of the noise amplitude A/" on the ratio r' and r. 

We see that if the amount of noise JV is high, it is possible that 
r' becomes negative (see Fig. 15.21b . which simply means that 
after factoring in the effect of J\f, quantum teleportation would 
perform worse than classical teleportation with r — and 
J\f — Q. More interesting is the case of broadband noise when 
JV — > JV{oj) contaminates the whole frequency range relevant 
for teleportation. In the same way that we had deduced Eqs. 
( 14.131 1 from Eqs. (|4.8l l using the wave-packet operator (14.1) , 
we define an effective noise level A/'eff — J fs {oj)Af{io)duj and 
write 

^out = ^in - V2er'''"vx - %/2 A/'eff y„, 
Pout = Pin + V2e-'''"wp + V2Ks zn, (5.20) 

so that finally it is possible to take into account the effect of 
Af{uj) by redefining rgff as in Eq. ( 15.211 ): 

J'eff ^ = reff - In yrrVV^J^ (5.21) 

In this section, we have shown how to estimate the effect 
of external classical noise sources on the output teleported 
modes. This simple model only works for Gaussian entan- 
glement and for sources of uncorrelated Gaussian noise. By 
Gaussian noise, we mean that the underlying quantum state 
used to trace out auxiliary quantum modes y and z is Gaus- 
sian. In this case, external classical noise becomes essentially 
indistinguishable from noise added by the teleportation itself 
due to finite squeezing. Furthermore our model is able to take 
into account any spectrum of noise JVito) by using the wave- 
packet mode function fs to estimate an effective level of noise 
A/^ff added to the teleportation. 

VI. BACK-TESTING AND CONCLUSION 

The first step to test the validity of our results is to check 
that the model of Sec. [Ill] we used for photon subtracted 
squeezed vacuum states works well with the experimental in- 
put states used in 1 1 1]. For that we need to estimate three pa- 
rameters: the squeezing parameter s, the loss parameter 1 — 77, 
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and the APD dark noise parameter e. In fTTl, a direct measure 
of the APD dark noise and event counts gives for e a value 
of 0.013. With the help of quantum tomography of the in- 
put squeezed vacuum state Ss\0) the squeezing parameter s 
is estimated to be 0.28. This tomography is done using the 
wave-packet function fs as a filter of the measured homodyne 
currents and without conditioning on the APD triggers. This 
means that s is actually an effective squeezing parameter in 
the sense of Eq. (14.12b . taking into account the bandwidth of 
the OPO used to generate the state S's|0). Finally, to estimate 
the value of rj, we use the equation ( 13.1 lb with the value of 
W°(0, 0) obtained from a quantum tomography of the input 
state and obtain 77 = 0.80. This value is slightly different 
from the one estimated in [11] due to the non-zero value of e. 
With these three parameters known we can numerically sim- 
ulate 14^'" using the results of Sec. |III]and compare it to the 
reconstructed Wigner function with the overlap formula, 

0{Wa,Wh) ^2TrJ J dxdpWa{x,p)Whix,p). (6.1) 

However, this formula does not work so straightforwardly in 
our case: if Wa and Wi, are mixed states and even though 
Wa — Wb, the overlap given by formula ( 16. lb will not be 1 but 
rather the purity of Wa- We therefore use a modified version 
of the above formula with a renormalization factor taking into 
account the purity of both quantum states: 



0\Wa,Wb) 



0{Wa,Wb) 



{0{Wa,Wa)0{Wb,Wb)) 



1/2- 



(6.2) 



With this modified overlap formula (16.2b we calculate an over- 
lap of 0.987 between our model and the reconstructed state. 
The L2 Euclidian distance d{Wa, Wb) defined by 



diWa,Wb) 



J J dxdp \Wa{x,p) 



Wbix,p)\- 



1/2 



(6.3) 

between the two states is found to be 0.05. Finally one could 
also choose to maximize O' rather than fitting the value of 
W^™(0,0) to estimate 77. However, because this approach 
can lead to a value of iy"(G, 0) significantly different from 
the experimentally measured value, we chose to directly fit 
W"(0, 0) instead. 

The second step is to estimate the broadband EPR parame- 
ter Teff, again choosing one of two possible methods. A first 
method would consist of directly measuring spectra of the 
EPR correlations between Alice and Bob and then using the 
mode function fs to obtain an estimation of r^ff. While this 
method automatically takes into account homodyne finite ef- 
ficiency and phase errors, it does not probe any imperfections 
of the classical channel. A second method would consist of es- 
timating r'eff with a measure of vacuum teleportation fidelity, 
with the added benefit of taking into account the whole pro- 
cess of teleportation. In the simple case of the vacuum state 
|G) as an input state, the teleportation fidelity and the EPR 
correlations parameter r are directly related by the relation 



To measure the fidelity i^teie we first make a tomographic re- 
construction of the teleported vacuum state. To specifically 
estimate r^ff we consider the wave-packet vacuum state |0) 
defined by 



|0)/. 



dLofsiio)\0), 



(6.5) 



1/(1. 



(6.4) 



As before this is simply done by using the wave-packet mode 
function fs as a filtering function in quantum tomography. 
With this second method, r^ff is estimated to be 0.795 in lUlll . 

Finally, the third and final step consists of checking our 
prediction of T4^™' and particularly T4^™'(0, 0) using all the 
known parameters. This is done by evaluating Eq. (13.12b and 
we obtain the value T4^°'-"(G, 0) = -0.0243, in agreement with 
the measured output negativity of —0.022 in [11]. As before 
since e 7^ the value of W^°"'(0, 0) estimated here is slightly 
different from the estimation found in [TT]. We also check the 
overlap O' between the measured state and the predicted state 
and find a value of 0.988. The L2 distance between the two 
states is 0.04. 

Overall, the results of Sec. |lll] are numerically in good 
agreement with the experimental results of ifTTll . Our model 
uses a set of only three parameters and even though more com- 
plex models for photon subtracted squeezed vacuum states 
exist, we found it was not necessary to use them. This con- 
firms the validity of our initial assumption to only consider the 
small R, small s regime. Our approach to multimode telepor- 
tation in Sec. |IV]has the double benefit of intuitively picturing 
broadband operations in term of wave packets, while at the 
same time allowing us to use the usual results of single-mode 
teleportation with simple renormalization prescriptions . 

In conclusion, we have developed an efficient yet simple 
framework to model the properties of multimode continuous 
variable teleportation. Although this work is tied to a specific 
class of non-Gaussian states, it is natural to ask if a similar ap- 
proach can handle more general non-Gaussian states. As it is 
known that any Wigner function can be approximated by suc- 
cessive displacements and photon subtractions, it is in princi- 
ple possible to describe any non-Gaussian states in a system- 
atic way that would be compatible with the Gaussian convolu- 
tions needed in Sec. |III1 A second and harder issue would be 
to identify criteria more robust than negativity to decide on the 
success of teleportation and Gaussian operations in general in 
the context of non-Gaussian non-classical input states. 
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